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Historical  Synopsiso  After  the  writer  develrped  his  theory  on  the  plastic  buck-* 
ling  of  plates  and  shells  he  applied  it  to  cases  of  homogeneous  stress 

distribution  for  the  most  important  cases  of  loading  and  boundary  conditions 
and  showed  its  excellent  agreement  with  test  results  fs],  fsj,  He  then 

Obtained  a  contract  from  the  Office  of  Ordnance  Research,  U.  S.  Army,  for 
applying  it  to  several  cases  of  non-homogeneous  stress  distribution  in  plates. 
The  results  of  this  work  was  published  in  three  papers  ^ [^]* 

In  an  earlier  paper  }  he  had  applied  his  theory  to  the  plastic  buckling 
of  a  cylindrical  shell  under  axial  compressive  load.  As  also  explained  in 
Reference  k  such  a  cylindrical  shell,  if  the  tangent  modulus  E  is  zero,  has 
the  same  plastic  buckling  stress  as  a  wide  column  with  the  same  half  wave 
length  of  buckles.  In  the  case  of  bending  of  such  a  shell,  the  plastic  fibers 
are  supported  laterally  by  the  elastic  one^  so  that  for  low  values  of  E^  the 
edge  buckling  stress  can  be  expected  to  be  much  hi^er  than  for  pure  compression 
vith  the  same  tangent  modulus  E  and  secant  modulus  E  .  Therefore,  the  present 
contract  was  applied  for,  in  order  to  compute  plastic  buckling  stresses  for 
several  cases  of  plastic  buckling  of  cylindrical  shelly  sub  jected  to  bending 
and  subsequently  derive  from  the  results  of  a  simple  design  formula  for  comput- 
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ing  the  plastic  reduction  coefficient  to  be  applied  to  the  elastic  buckling 
stress. 

Method  Bmployed.  Prom  Reference  3  the  differential  equations  for  the 

plastic  buckling  of  a  cylindrical  shell  under  axial 


compressive  stresses  are 
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where  u^,  v,  and  w  are  the  displacements  in  the  axial  (x-),  tangential  (aO-)^ 
and  radial  (z-^  directions,  respectively,  a  and  t  are  the  shell  radius  and 
thickness,  and  A,  B,  D,  and  F  are  plastic  reduction  coefficients,  which  are, 
for  example,  given  in  reference  k  as  functions  of  the  tangent  and  secant 
moduli  and  of  Poisson's  ratio  in  the  elastic  range.  These  equations  may  be 
reduced  to  a  partial  differential  equation  in  w  only,  namely 
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where 
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Eq*  (^)  was  reduced  to  an  ordinary  differential  equation  by  the  substitution 

where  Y  is  a  function  of  y  =  aO  only  and  ^  is  the  half  wave  length  of  huckles 
in  the  axial  direction.  Dividing  the  circumference  of  the  shell  into  32 
spacings  this  ordinary  differential  equation  was  written  as  a  S6t  of  17  (due 
to  symmetry)  finite  difference  equations ,  using  second  order  finite  differences, 
as  were  also  used  for  the  case  of  buckling  of  a  plate  under  plastic  bending 
in  its  plane  (reference  5)-  This  results  in  the  following  difference  equation 
for  an  arbitrary  point  ks 
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In  here  h  is  the  spacing  ita/l6  of  the  pivotal  points  on  the  circumference^ 

A  -  where  ^  is  the  half  wave  length  of  buckles  in  the  axial  (x- ) 

direction.  Further  =12  0^/Et  where  0*^  is  a  reference  stress,  and 
^  •  Th^  17  equations  (5)  lead  to  a  17th  order  determinant  which 

has  to  be  equal  to  zero  from  which  the  edge  buckling  stress  can  be  calculated. 

Cases  1  and  2  were  calculated,  where  the  plastic  range  extends  to 
(5/l6)jt  and  (3/l6)7t  from  top  and  bottom,  respectively.  For  the  plastic  range 
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several  constant  values  of  the  tangent  modulus  were  ass\iraed  and  also,  to 
establish  the  optimim  half  wave  length  £ ^  several  ratios  ^/a.  The  diffi¬ 
culties  with  this  project  were  (l)  to  find  graduate  students  that  were 
interested  in  computing  the  17th  order  determinants  that,  at  the  same  time, 
worked  sufficiently  accurate,  (2)  the  fact  that  silll  inaccToracies  occurred, 
which,  in  most  cases,  leads  to  non-converging  solutions  of  the  determinant, 

(3)  the  students  cannot  spend  very  much  time  on  this  work,  so  that  it  progresses 
slowly  and,  after  the  student  graduates,  a  new  one  has  to  be  found,  (4)  The 
difficulties  of  finding  the  numerical  errors,  and  (5)  the  evaluation  of  the 

determinants  for  finding  the  lovrest  eigenvalue. 

\ 

Only  the  last  difficulty  was  finally  solved  by  having  the  determinants 
evaluated  at  Bell  Aircraft  Company  in  Buffalo,  New  York,  on  their  704  IBM 
machine.  The  result  was  that  only  the  following  cases  could  be  evaluated, 
where  the  last  column  givej  the  edge  stresses  0*^  at  incipient  buckling.  It 
was  assumed  that  a/t  =  1000/  . 
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Table  I 


Case 

<^/E 

1 

.5 

.00000978 

1 

.7 

.0000505 

2 

.5 

.0000373 

2 

.7 

.0001137 

1 

.5 

.00206 

1 

.7 

.001905 

1 

.9 

.001984 

2 

.5 

.002093 

2 

.7 

.001944 

2 

.9 

.002061 

.7 

.002032 

1 
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In  order  to  establish  a  general  design  formula  also  values  for  E^/E»^ 
-^/a  =  .3  and  several  cases  for  E^/E  of  .1  and  .5  have  to  be  evaluated. 
Determinants  for  most  of  these  cases  were  computed  (some  after  the  contract 
was  finished  at  my  own  e:qpense)  but  they  will  still  have  to  be  checked  for 
accuracy.  This  will  be  done  if  time  for  it  is  avialab le. 

The  numerical  results  obtained,  as  given  in  Table  I,  sure  not  sufficient 
for  deriving  conclusions  for  other  cases. 
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